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ABSTRACT 



This study presents the results of an experiment in objective 
analysis of oceanographic data for a limited area. The objective 
analysis is designed to provide a reliable operational system for 
tactical use in coastal waters. It is shown that this approach makes 
it possible to obtain a very detailed analysis with good vertical 
consistency and that only a relatively small amount of highly accurate 
data is required. The procedure includes a polynomial interpolation 
and a dynamical interpolation. The dynamical interpolation which is 
introduced involves the analytical solution of a finite differential 
equation by means of a rotational computing molecule. The procedure 
requires only a small computer and little computer time. This method 
will provide a basis for short-time forecasts of oceanographic para- 
meters using only small computer centers or even time sharing systems. 
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I. INTRODUCTION 



Analysis of oceanographic data has been traditionally a very slow 
and essentially a subjective process due to the time involved in an 
oceanographic expedition and the small amount of data obtained. In- 
creasing the number of systematic observations in time and space 
requires faster analysis procedures and objective methods. Such methods 
will provide information for regions of sparce data. The analysis pro- 
cedure should also standarize sequential analyses so that the results 
may be* compared. 

Objective analysis techniques are currently used on an operational 
basis by the Navy and other government agencies concerned with atmos- 
pheric phenomena. Background information available for formulating 
an objective analysis scheme for oceanographic data is primarily that 
from investigations and from experience with schemes for meteorological 
data. However, there are significant differences between both the 
type and distribution of meteorological and oceanographic data which 
have to be considered. 

The primary purpose of most objective analysis schemes is to 
transform data from observations at irregularly spaced points into data 
at points within a regularly arranged grid (grid points). A subjective 
method for doing this would be to have an analyst manually (subjectively) 
draw isolines and then interpolate (objectively/subjectively) to obtain 
values at the appropriate grid points. An example of an operational 
grid is given in Figure 1. 
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FIGURE I. • Example of the grid used by the Joint Numberical Weather 
Prediction Unit (After Cressman 1959) 
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Procedures used in objective analysis have to be systematic and 
quantitative so that they can be performed by numerical computations 
on a computer. Computer processing is by its very nature, objective. 

In practice, once the observations are relayed from reporting stations 
all procedures leading to an analyzed chart should be completed within 
a minimum of manual intervention. 

A. BACKGROUND 

Research and the subsequent development of objective analysis 
methods took place in response to the need to prepare large numbers 
of meteorological observations representing the atmosphere, both in 
the horizontal and vertical, for numerical weather prediction procedures. 
Numerical weather prediction requires, in addition to values at reg- 
ularly spaced grid points, editing (removal of erroneous observations) 
of these data and checking for consistency and continuity in time and 
space. These requirements would also apply to objective analysis 
schemes for oceanographic data. 

Once a grid of equally spaced values, checked for possible errors, 
is generated by an appropriate objective analysis scheme further work 
can be performed readily on a computer. These steps also could be 
performed in subjective analysis but the time required to obtain the 
completeness achieved by an objective analysis scheme would be out of 
proportion to the time needed for the computer analysis. 

An early publication on the subject of objective analysis of geo- 
physical data was by Panofsky (1949). Various approaches for per- 
forming objective analysis with meteorological data were presented 
during successive years by Gilchrist (1954), Bergthors on and Doos 
(1955), Johnson (1957), Cressman (1959), and Smith (19o2). The latter 
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references describe versions of operational objective analysis schemes 
presently used in the principal weather centers of the world. However, 
most of the currently used objective analysis schemes are based on the 
"Operational objective analysis system" described by Cressman (1959). 

The general approach in these schemes is a very straightforward 
corrective approach in which an initial (preliminary) field is modified 
by available observations and the result is then used as preliminary 
field which is again modified with the observations but with more 
resolution. The repeating corrections are called scans. This process 
may be repeated many times but in general satisfactory results are 

obtained with the first or second scan. 

In the corrective approach, an analytical model defines the cn- 



terian for checking the validity of observed values and determining 
values at grid points. Determination of grid point values from obser- 
vations at nearby 1 stations is achieved by some form of interpolation. 
One method -of interpolating for a scalar field is to represent the 
contours of isolines surrounding the grid point as a function of 
position. This may be accomplished by means of a surface least squares 
polynomial interpolation. The order or degree of the polynomial in the 
least square interpolation depends on how detailed the final field 

should be; 

"The degree of smoothing will depend upon the nature of the fore- 
cast and the particular model selected Smoothing to eliminate 

unwanted, though real, wave lengths has been taken to fall within 
the province of the forecast rather than that of the analysis. 
(Johnson 1957) 



1 Nearby is normally considered to be within one grid distance. 
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Oceanographic analyses have been made only as a related factor in 
weather analysis. Climatological forecasts of ocean surface temperature 
anomalies have been made by Adem (1969), Jacob (1967), and Adem (1970) 
using the equation for the conservation of thermal energy. 

The type of data expected in an oceanographic survey introduces 
the problem of holes between transects and therefore the need of some 
type of dynamical interpolation. Typical cruise tracks appear in 
Figures 2 and 3. Provisions for lack of data for meteorological 
analysis of a given region when good data is available at the boundary 
have been studied by Thompson (1961), Richardson (1961), Smith (1961), 
and Smith (1962). The method generally used is to "advect good data 
from outside the region (or "hole")" (Smich 1962). 

S. PURPOSE AND SCOPE OF THIS STUDY 

The purpose of this study is to develop an objective analysis 
scheme for oceanographic data which will meet operational requirements 

described in the following paragraphs. 

The high cost or operating oceanographic ships at sea is a primary 
reason to develop methods to complete instant detailed analysis. Such 
an analysis would show interesting features based on early observations 
but not considered during the initial planning which would permit 
scheduling of complementary observations during the same operation at 
sea. Furthermore, operation requirements of submarine and antisubmarine 
warfare have shown the need of faster and more accurate short term 
forecasts of oceanographic parameters which could be produced through 
the use of numerical procedures. 

The scope is to develop a scheme which would resol'-e fine details. 
Small thermal inhomogeneities have great effect in sound propagation 
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FIGURE 3. Cruise track showing station locations along Region II of Caribbean coast of Colombi 
(After Sanchez 1969) 



used. Those fine 



especially when relatively high frequencies are 
details cannot be observed nor forecasted with large scale schemes and 
therefore, this study considers small scale objective analysis methods 
for oceanographic data, the method will include some type of dynamical 
interpolation. The scheme will be evaluated with respect to its use 
for tactical purposes and as a "feed back" for planning during ocean- 

ographic surveys. 
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II. approach 



Data obtained from the standard oceanographic stations differ from 
the meteorological data in the number of stations, the accuracy of 
observations and in the fact that the former nearly always has to be 
considered synoptic 1 for long time steps over small areas. In general, 
an initial field will not be available in oceanographic analysis nor 
will there be enough data to obtain a complete interpolated field. 
Therefore , new procedures have to be formulated in the development of 
the objective analysis for oceanographic data. Techniques described 
above for meteorological analysis will be applicable only in a relatively 

small part in this particular problem. 

The objective analysis developed is a numerical scheme in which a 
least squares surface polynomial is used for interpolating and smoothing 
the available data to nearby grid points. Then a dynamical inter- 
polation is used to compute values at those grid points which are in 
"holes" or regions of no data. The dynamical interpolation is achieved 

by means of a rotational molecule. 



A. POLYNOMIAL INTERPOLATION 

The first step in the scheme is to reduce the original data to 
neighboring 2 grid points using surface least squares polynomial 



1 Synoptic data represent observations taken simultaneously 

2 Neighboring refers to data within one grid distance. 
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interpolations . A second degree polynomial is used when sir. or more 

data points are available in a given square around the grid point, and 

fits t degree polynomial when only four or five points are available. 

A weighted gradient propagation is used when only two or three data 

points are available. When only one observation is available in the 

entire square around the grid point some criteria about the distance 

is established and if this distance is, say, less than five mile 

•n h P t-o give the same value to the grid point, 
best approximation will be to give 

otherwise this value is neglected. 

The polynomial solved for the second degree polynomial interpolation 



is : 



♦ a to x + a„ 




+ a ie r+ajj+a 



CD 



The 



first degree interpolation formula is: 

x + a„xy + 



( 2 ) 



and the weighted gradient propagation 



formula is: 



z.-z*w 



Cx,-x. )(y,-x) 



( 3 ) 



Where * 4.0 ° £ ““ *** P ° iat C ° nSidered - 
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The matrix obtained in the first two cases is solved for 6[ 0o , 

since the grid point is at location 0,0 , using the subroutine SIHQ 
from the computer library of the Naval Postgraduate School. 



This scheme for obtaining values for some grid points enables the 
maximum use of data. 

In meteorological analysis the regression has been used and also 
the correction method. In the latter, it is assumed that the difference 
between the preliminary value and the observed value at a given station 
is equal to the difference between the derived value and the preliminary 
value at the grid point. Another approach in the correction method is 
to assume that the gradient of the preliminary field at the grid point 
is representative for all the area between grid points and the station 
(gradient propagation). Then the correction is a function of the 
distance between grid points and station. In these procedures a 
statistically determined weight function (weighted gradient propagation) 
is required. 

In all interpolation methods, the distance between observations and 



a grid point becomes a crucial factor; for example, 

"Quite often it is not possible to get a reasonable analysis only 
by means of interpolation between synoptic observations. It is 
quite clear that the distance between observations must be sma 
compared with the size of the system to be analyzed." 
(Bergthorsson and Doos 1955) 

Therefore, distance limitations introduce the need for a dynamical 
interpolation whenever "holes" are present. 



B. DYNAMICAL INTERPOLATION 

The initial step in developing a suitable dynamic interpolation is 
simplification of the basic conservation equations for a scalar quantity 
in a fluid medium. In order to apply a simplified equations to a set ol 
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regular horizontally and vertically distributed data, certain assump- 
tions are required about the flow characteristics. These assumptions 
have to be based on conditions observed or generally accepted for an 
oceanographic environment. The simplification of these equations, the 
assumptions involved, and the numerical procedures used to apply the 
final form of the equation are discussed in the following paragraphs. 

1. Analytical Model 

Consider the analysis of a conservative scalar field. The 
available data will represent quasi -synoptic observations provided by 
standard stations in a relatively small area of the ocean. Therefore, 
some of the following information could be available for each station: 

(1) Temperature 

(2) Salinity 

(3) Wind velocity 

(4) Dynamic height anomaly 

(i) Sigma-t 

(6) Sound velocity 

at observation depths and at standard depths. 

An analytical model for analysis of temperature, salinity and 
sound velocity would be based on the following steady state conservation 
equation for a non divergent fluid: 



U '& Dy 2AD-A & 



K 



m 

2z. 



(*) 
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where 



u, v, and w are the components of the velocity vector 
K is the eddy viscosity coefficient 

Studies have shown that the following approximation is possible 




where; 

"The condition for neglecting the horizontal terms is that they 
approximately cancel each other. Since those terms have been 
shown to be of the same order of magnitude and opposite signs, 
this condition is possible." (Overstreet and Rattray 1969) 

Assuming condition (5) the conservation equation becomes: 




+ Of' 



d_D - d 



^ r 



" dx 






"h r 

+ 4-\ k. 



si{ 






( 6 ) 



Assuming isotropy and horizontal homogeneity the two dimensional 
Equation (6) becomes: 



7 / 2 O + IT 20 - jc [20 + &P 



(7) 



or 



ji 

K„ 
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where 



P - 

u 

P can be computed for each grid point if the direction of the current 
at each grid point is known. However, unless direct measurements of 
the current have been made near each grid point this is not generally 
the case. Therefore, the next best solution is to compute the current 
for each station using available data and extrapolate or interpolate 
it to the nearest grid points assuming that the current is uniform m 

the area. 

The current consists of four components 

V- Vj + Vj+V^V, (io) 

where 

Yj is the drift current due to wind action 

Y^ is ’the geos trophic current computed from the density field 
Vi is the baro tropic current 

K is the tide current 

If a strong baro tropic current, Vi , i- s present m the region to be 
analyzed, climatological data would have to be considered for its 
definition. 

Tidal currents may be significant because oceanographic data 
is only quasi-synoptic as noted before. However, currents considered 
are generally time averaged for the entire time of all the observations 
Therefore, observations would represent a time lapse of two or three 
days or at least 24 hours. It is possible, therefore, to neglect tidal 
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_ However, when the observations are truly 
currents in the average. However, 

synoptic, tidal currents must be considered, i. e. 

"Tidal currents are simply are peri- 

forces . Since tidal f^^y 

odical and these current 8^ considered as part of the ad- 

water , tidal currents will^^ st „ cture Forecastrng, ASWEPS 
vection problem. 1 UCI e ° l 

manual series, v. 5, p. *■*■) 

Hereafter we neglect tidal and barotropic currents 

horizontal components of the current vector become. 

U " Mel + 



qr ■= ITd + V) 



(u) 



and 



p _ qjjutM 

r & * 



Similar formulae have been suggested by Arthur (1966) and Aden, 

(1970). _ . , 

These velocity components can be computed from the data avail- 
able in a standard oceanographic station. 

Considering the drift current, we examine possible relation 

between it and the obsereved wind, benoting the surface wind velocrty 
as Va, the following formula apply for Va along the 



_^7f/o)2 7T -A 

io-ve cm ( 4 s - o i ) 



(12) 



-(%) z ,. 2z ^ 

^ s Vf 0*) 
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where 



V = 0.0 253 



tim 



and 



J) _ 3.6? (ML 



V = JULLZJkl^ 
' " jpUrtf 



S 



or 



IV I- 0> ort/ttet. 



D = 



iv, I 



7.6 



f 



or 



|V.| > I i rm/jM. 



(P 1 



is the latitude 



For the general case, when 



the Y-axis is not along the wind direccici 



(13) 



Equations (12) becomes (Adem 1970), 

If - Qs o } AA f V.* 0 - V» Am O') 

* \CAm~(f V 

/y _ ojiiiA ( OJl 0 - V a Aim & } 
d jfJAAf ' 

Where $ is the angle between the surface current and the surface 
wind. 0 is positive when the current is to the right of the 

surface wind. 

The geostrophic current ( ^ )> between two stations 

may be computed from the dynamic height anomaly using the Ekman 
approximation: 
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f K - f« ’ f M 

.# -if. . ® %? 



( 14 ) 



where / is the Coriolis parameter 

is the geopotential anomaly 

Assuming that the density increases steadily with depth, the relative 
geos trophic current between two stations A and B, can be expressed as 



/ir- % -- ^(aD a -^A) 



(15) 



where L is the distance between A and B 

A further simplification of Equation (11) may be achieved 

assuming that the geos trophic current has no contribution to the 



advection of density anamolies; 

"The major contribution to advection (of heat) comes from the 
wind drift current since the slop current is assumed to be 
geos trophic, or flowing rarallel to the isotherms. That is, 
there is no advection (of heat) by water flowing parallel to the 
isotherms; only cross section flow causes a change m the 
perture at a point." (Ocean Thermal Structure Forecast, SP-105, 
ASWEPS manual series, v. 5, p. 33) 

If this approximation can be accepted, Equation (11) becomes, 



/y - /vd 

(16) 

U-- Uj 
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p = 



a C l 

U, 4 



Therefore, many degrees of approximation may be used to compute 

P and they can be summarized in three cases: 

a. Currents at the station are known . If the true current has 

been measured at each station, P may be computed directly or if the 
current has been computed, P may be computed from Equation (11). The 
same P may be used to nearby grid points provided the changes in current 

direction are not too large relative to the grid scale. 

b. The wind velocity is known but t h e true current iS i not 
known for each station . In this case only the drift current can be 
considered and Equation (16) must be used to obtain the direction of the 
current for the grid points using the wind direction nearby. Also, 

in this case empirical consideration may be used for allowing to. tne 
drift current. These empirical considerations are: When slope and 

gradient currents (due to mass transport by the wind) have been set up, 
they affect the direction of the original wind driven current, and when 
a steady state has been reached, the deflection of the surface from 
the down wind direction will be reduced to near 20 degrees at low 
wind speeds or to less than 12 degrees at high wind speeds (Summary of 
Technical Reports of the National Research Committee, v. 6A) , (Ocean 
Thermal Structure Forecasting, SP-105, ASWEPS manual series, v. 5), 
(Hubert 1964). 

C. Only the wind field is known but oth er factor may be 
important . This situation may arise when the coast or bottom topography 
would invalidate the assumption that the current is at a constant angle 
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to the wind direction. For this case it is possible to use the follow 
ing relation to obtain as estimate of the direction of the current at 

each grid point « 



V 



(Kitaigorodskii 1961) 



and therefore, 



m - COS TAN T = RATIO 

K 



where V is the fluctuating velocity vector. 

However, sufficient surrounding data are required to apply a numerical 
solution for obtaining these variables. In order to compute ±VL 
a P value, from the wind field is assumed for each grid point as m 
b. or for the whole grid if the wind field is more or less uniform. 
Using the assumed P, values of U/K are computed for all possible grid 
points using Equation (8) . The numerical sheme to perform this will 
be explained in the next section. 



Defining 



\ L - U. 
k " K 






and assuming that eddies present in the area may introduce a great 
range of directions in the fluctuating velocity and given that 1^- 
is a constant, the maximum value of U/K should correspond to a minimum 
value of • Also, if ^ is small enough, it is possible the 

approximation 



I V] -u Li maximum “ RATIO 

K K 
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This approximation can ba made if it is assumed that the mean current 
is thermohaline and therefore flows parallel to isotherms so it does 
not contribution in the advection of heat.i. e. 

" + V' VT = V*T 



^ vf‘ 



if V ±VT 

Therefore, only the fluctuating component is obtained when 
and VT are computed from a field of temperature to solve for 

V . For each grid point at which U/K is computed a/tf/K can be 

obtained using the formula 




and the new P 



V/X. - or 

U/k u 



Using this U/K and the new P Equation (8) can be solved for T at next 
grid point if enough nearby grid values exist to use the rotational 
computing molecule presented in the numerical model (seccion 2). 

The procedure above can be repeated many times but it has been 
observed that usually the first approximation is satisfactory. Because 
the quadratic equation (17) has two solutions, the wind field is used 
in the selection of the correct solution. This choice can be objective 
and included in the computer program. 
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2 . The Numerical Model 



Knowing P, U/K can be computed whenever the parameter is 
known in four neighboring grid points and then, using this value 
U/K and the corresponding P the next unknown grid point parameter can 
be computed. The problem reduces to solving the differential agnations 



u 



'H L 



20 +v 



= Id 

2 x z 






( 8 ) 



uj2u. - IlD 

T "32 32* 



The corresponding Finite Differential Equations may be solved using 
the same scheme. 

■ The following is a description of a direct solution of the 
partial differential equations by means of a rotating computational 
molecule in such a way that a solution is possible oven though data are 

scarce . 

The horizontal continuity equation, Equation (8), may be 
approximated by the following computational molecules: 



(1) a. 



Z f (XO = -2 (K,L*2.)-Z(K« t L)-ZmMZ) 



(i pr} 

K'XR 



fa 



Z(KH,L*1)~Z(K*B,U2') 
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b. 



t )/ 2 ■ -j 

7 ( KL \- 7 \ Z ( XM 2 )- Z («+ 2 , L )\ 

( ' ^ 2.0*K*~XR L J 




+ 2.0*2 (K^O-zCKs-tz)- p^~ h(xMO-Z(KH$l 

K * ysc [ j 






( 3 ) a. Z p ( K j L ) - ^.0*2CK-1 > L-tl)-2(X t Li2')-2(K-2 J i)-Z{K-2 J Lt2) 

_ fU ’ iV ’ f 1 hex -/, up - z ( x-z *■**)! 

L J 



b. 



z w = Z/*- 1 ) - 



(U'tnr r f l \jukLfz-\ - 2 Oz,/) 

2 .o*k*"ar L 



( 4 ) a. Jf - " (T^^]zLK-l,L + l)- 2 . 0 *Z(K- 1 t L)iZ(X- 1 A-t^ 

+ 20*2 (K- 1, L) ~ 2 CK- 2 , L)- JL- \z(K‘t,L)-' 2 (Xr 2 J L) 

K*xscL 

b. ZCW ^ ftO - 



( 5 ) a. 



2 ^ = 4 o*z(K-i M - , )- z( *'*' L )~ z(Kt ~ z) ~ z( *' 1 ' L '* ) 

_ (u 2 i/r*) H \z(K-i, i-n-z 

1 J 
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b. 



(6) a. 



b. 

(7) a. 



b. 

(8) a. 



b 



(9) a 



ZM- Z,(v)-$£jg£ [zrv- 0 -zr*i,o] 



2 r (fify = -(gj[z(*'i l-<)-2-0'Z (V-0+ 2 f/t-Ji-J j 
+ 2.<?» - Z(Ki^)-£~Jz(Ki-0-Z(K^)j 



U 



7 (k L ) - 7 ( k t i)- 

k- K ‘ J £? K ' ' 2.0*K+XSC 






J p (X t L) = 40*Z(Kit J L-l)-Z(Ki2 > L)-Z(K J L-2)-ZO«27-2) 

s</2 r 



_ X /&0 

K-XR 



Z(K>l,L-0- Z(KH,t-t) 

' M Z(k*',l)-Z (K,L-Q 



7 (k 0 = 

2 r ~] 

~7 (X t L)- ^^2CKtl,L*i)-2.D*Z(X+ij L ) + Z(K*2jL ) j 

4 2 . * Z(.Ki1 J L)-2(K+2 J L)- JA— f Z (KHjO-ZO**! 1 )] 

K*X£C L J 



7 f*/) r 7 (V) 

2.G*K» 



2,c*x»y£C 



ZC^t.L^O'ZikH^-i) 



' U^T^J 



ZftO = Z r (X,l)-[z(*>tO-Z<*-W] ]ZAxsc + 



2.0» 
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in each case the competing molecule has been divided in two parts in 
order to give more flexibility to the computation and be able to compute 
the Lapacean and the gradient separately. In order to be able to use 
all molecules, the X-axis scale, XSC, must be equal to the Y-axis 
scale, YSC. If the proceeding is not the case, local distortion may be 
introduced. Also the above procedures are not iterative and, therefore, 



stability problem will not arise. 
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FIRGURE 4. Sequence of rotation of the computing mole ale around 
a grid point with no value U(K,L). 



32 




■i — q —a —a — 

2 

^ ••••#•• 
4 — 0 — 0 — 0 — 0 — 

5 




Un known (K,L) (q 



HRE 5. A complete rotation of the computing molecule. 
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HI. DESCRIPTION OF DATA 

In order to verify the proposed scheme, data were selected from 
areas with good data coverage and well known physical features. 

One area considered was the California coast between San Francisco 
Bay and Monterey Bay and the corresponding data were obtained by the 
AGOR 1 cruises of September 1969 and 28 April - 4 May 1970. These data 
correspond to standard oceanographic stations (Figure 2). However, 
not all the information was readily available due to the fact that the 
data were reduced for specific uses. The data available were; 

a. Wind 

b. depth 

c . temperature 

d. Sound velocity 

e. Beam transmittance 

at observation depths only. Under these circumstances the general 
procedure. Equation (11), doesn't apply and therefore, the modified 
procedure, c. , was used. 

A second area was considered in order to verify the objective 
analysis procedures at low latitudes. This second set of data was 
obtained by the Colombian Navy's Division of oceanography in the 
Caribbean sea off the Colombian coast. In this case all the data from 



Cruises by ships of the U. S. Oceanographic Office under inves 
tigation plans by the Naval Postgraduate School. 
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a standard oceanographic station were available. However, the same 
simplified procedure, c., was used in order to obtain results com- 
parable to those for the California coasts case. 



35 



XV. DISCUSSION OF RESULTS 



A. COMPARISON OF OBJECTIVE AND SUBJECTIVE ANALYSIS 

A subjective analysis, Shepard (1970), of sea surface temperature 
for the area along the California coast appears in Figure 6. Shepard 
performed the analysis in order to make a comparison of oceanographic 
parameters during upwelling. The objective analysis of the same data 

appears in Figure 7 . 

For the objective analysis, Figure 7, the grid was rotated 47 degrees 
in order to follow the general trend of the coast. Because the mean 
wind was generally along the X-axis during the cruise, an initial 
value of P equal -0.27 or a current 20 degrees to the right of the 
X-axis was selected. 

A comparison clearly shows greater complexity and more details in 
the objective analysis than the subjective analysis. The obvious 
question that arises is; do these details represent artificial additions 
or real features with physical meaning? Also, do they represent features 
expected with the inherent assumption of the objective scheme? 

It may be stated that subjective analysis is a linear interpolation 
which normally depends only on the available data for one oceanographic 
variable at a time. This generally excludes consideration of currents. 
Representations in intermediate regions with no data are generally 
smoothed contours. The resulting scalar field (as depicted in Figure 6) 
is very familiar in the oceanographic work and a departure from this 
picture is often viewed with suspicion. Adverse appraisal of a more 
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FIGURE 6. Subjective analysis of sea surface temperature for 28 April - 4 May 1970 
(After Shepard 1970) 



\ 
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FIGURE 7. Objective analysis of sea surface temperature for April - May 1970, showing wind field. 
Arrows are the V/K = U/xi + J vector. 
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FIGURE 8. Comparison of subjective and objective analyses of sea surface temperature for 
California coast. 28 April-4 May 1970. 



detailed representation is not totally justified as stated in the 
following passage: 

■'Linearity between widely scattered data should not be assumed. 

It is a paradox of analysis that the less data the analyst has 

Che more the ^lysis show^th ^^Slnd P fhe presence 

of'dftaUed structures becomes apparent but not entirely 

by the data that the analysis becomes difficult. (Ocean Therm 
Structure Forecasting, ASWEPS manual series, v. 5, p. 23). 

Other features present in the objective and absent in the subjective 

analysis are the circular patterns or pools of relative warmer or cooler 

water. This type of temperature distribution is intuitively logical 

as most swimmers can tell and physically may be expected. Subjective 

analysis procedures are not suitable for detecting them because they 

do not normally involve considerations of temperature along with 

advection currents and diffusion terms. 

»t.jj .- pc have seldom been tuny delineated . - -~o 

of temUature observations crosses what may be a roughly or rcular 

Defense Research Committee, v. 6A, p. / ) 

The importance of including the influence of the wind action, i. e. 
wind currents, within the scalar analysis is emphasised further in the 
following physical arguments. 

Wind and wave action have a net result of mass transport of water 
and therefore both are capable of producing zones of convergence or 
divergence accompanied by a rising or lowering of the thermocline or 
may be responsible for water of one temperature being carried over 
water of another. (Summary of Technical Reprots of the National 

Defense Research Coiranittee, v. 6A) . For example, 

■'The effect of the wind in maintaining or deepening an isothermal 
surface layer is complicated by the fact that the wind sets up a 
current near the surface which may remove the mixed -ater from 
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, cl . fl _ fast as it is formed." (Summary of Technical 
Reports "^f ^national Defense Research Committee, v. 6A, p. 42) 

B. interpretation of results in terms of the analytical model. 

Additional interpretations of the features within the objective 
analysis are possible by considering the basic processes included in 
the analytical model. The physical description of the analytical model. 
Equation (4), used in the objective analysis scheme is that when a steady 
state is reached the advection of the scalar property is equal to the 
eddy diffusion of the scalar in the medium. With this physical inter- 
pretation, we can explain certain features of the circular distribution 
of properties when these features are frozen in time (stationary). 

These features are: 

1. Circular patterns, 

2. Tongue type patterns, 

3. Gradients 

4. Siz-e of closed isolines. 

With respect to the circular patterns, heating and cooling of the 
sea surface is not uniform for the entire layer of the analysis and 
therefore in any event a local circulation will be generated in order 
to obtain the equilibrium of mass. Warm water will flow from the 
source toward the periphery in all directions and the Coriolis accel- 
eration will deflect the paths to the right (in the Northern Hemisphere) 
until they become circular. This situation will exist as long as a 
disequilibrium in. mass exists and will become a steady state when 
advection out of the circle is equal to the inflow of warm water, a 
similar situation is reached in a cold pool where warm water from the 
periphery is going to the center of the pool following a circular path. 
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Decaying drift currents will also produce the same effect because once 
the wind ceases the drift current continues to move for a certain time 
under its inertia while being acted upon by the Coriolis acceleration 

and while being damped by the diffussion term. 

With respect to tongue type features, the oceanographic literature 

has abundant examples of tongue-like diffusion of properties such as. 
salinity and temperature. For example, it has been observed that in 
areas of convergence, convection cells, the surface water moves down- 
wind at a considerably more rapid rate than in areas of divergence and 
therefore the flow of water in response to the wind should be uneven. 

With respect to the gradients, because vertical mass transport 
occurs in the convective layer, it is clear that the assumption of a 
horizontally homogeneous mixed layer above the thermoclme is an over 
simplified concept and relatively high gradients may be possible in- 



short horizontal distances. This is especially true in the area 
examined which is in coastal waters, where an irregular wind field 
may have existed and where upwelling is known to exist. The analytical 
model was, of course, a non divergent and this may not have been the 
case. Therefore, greater than normal gradients inside the inertial 
circles may reflect over compensation for the upwelling or downwelling 
effect in the temperature. Although, this is, perhaps, an unrealistic 
result, it could be smoothed out by computing the amount of this effect 
using the grid points over nearby coastal regions (0. S. Fleet Numerical 
Weather Facility. Technical memo No. 5, January 1965). On the other 
hand, artificial intensification of the gradient may still be used to 
identify convergence and divergence zones by allowing a maximum gradient 
in a given region depending of local conditions. 
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With respect to the size of closed Isolines, the radius of the 
inertial circle for the type of circulation mentioned above varies 
with the velocity and the sine of the latitude: For a current of 

average velocity located in mid-latitudes this radius is about 15 miles 
which can be greater if the circulation is influenced by a horizontal 
component of gravitational force or other lateral force (Summary of 
Technical Eeprots of the National Defense Research Committee, V. 6A) . 
It is interesting to note that 15 miles is the order of magnitude 
observed in the circular features of the objective analysis under 



consideration . 

Therefore , it seems possible, with this analytical model, to relate 
circular and elliptical features to inertial currents due to differ- 
ential heating and cooling, convergences, divergences, decaying wind 
drift currents; tongue like features to normal diffusion processes; 
and intensified gradients inside circular patterns to strong upwelling 
or downwelting. However, gradients are represented very well within 
the analytical model. Reasons for the latter feature are however, known 



C. COMPARISON OF OBJECTIVE ANALYSIS RESULTS FOR DIFFERENT LEVELS 

In order to verify vertical consistency, a mixed layer of at least 
10 meters was assumed and an objective analysis of the 10 meter 
was made and is shown in Figure 9. However, the wind field was not 
uniform and, there was not enongh information about fetch and duration 
with available wind observations. Therefore, we do not expect to obtain 
an exact correspondence of the features of the two layers over the 
entire grid. Nevertheless, vertical continuity is observed when the 
main features of the two analyses are superimposed in igure 10. 
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Maximum and minimum temperature regions are marked with letters A 
through K in Figure 7 and apparent corresponding features in Figure 9 
have been marked with same letters with a prime sign. The same regions 
are shown in Figure 10 but delineated by only one or two isotherms in 

order to simplify the resulting picture. 

The comparison of thermal structure at the two levels shows an 

almost identical correspondence between regions 

A and A 1 
B and B f 
D and D 1 
E and E 1 
F and F 1 
H and H' 

I and I 1 
T .T 1 



This may be expected from good vertical mixing and more or less constant 

advection in the entire layer. 

Regions C and C 

K and K' 

G and G f 

have horizontal relative displacement from 5 to 12 miles. This may be 
due to a different velocity field for the two levels. Past history 
(i. e. two or three weeks earlier) would show how these advective 
distortions were generated. Similar results may be obtained for 

September data (Figures 11 and 12). 

Considering that each analysis is completely independent from the 

other, this comparison shows convincing verification c verti 
consistence . 
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28 April - 5 May 1970. 



D. OBJECTIVE ANALYSES OF SURFACE TEMPERATURE, SALINITY AND SIGMA-T 

FOR LOW LATITUDES 

Objective analyses of temperature, salinity and sigma-t for a region 
with mean latitude of 10° north appear in Figures 14, 15, and 16. The 
main reason for considering these analyses for a region not very well 
studied from the oceanographic point of view is to compare the effect 
of latitude in the size of those features interpreted to be inertial 

circles in the preceding discussions. 

The region of California coast considered has a mean latitude of 
37° north and feature A in Figure 7 appears to be representative one for 
this region. It is observed that its mean radius is about 8 nautical 
miles. The radius of inertial circles is inversely proportional to the 
sine of the latitude and therefore it would be expected that the radius 
of the inertial circle be increased by a factor of approximately 3.4 
when translated from latitude 37°N. to loV a corresponding feature 
in Figure 14 was marked A" has mean radius of about 35 nautical miles. 
Therefore, the test for latitude effect seems satisfactory from this 
gross comparison. However, more precise computations would be required 

before a conclusive result may be obtained. 

The region off the Colombian coast, considered in the low latitude 
study (Figure 13), receives the run-off of three main rivers: the 

Atrato in the south-west region of the grid and the Smu and Magdalena 
rivers in the east. Features L,M and N in Figure 15 correspond to places 
where the run off occurs. The patterns observed at these points coin- 
cide very well with what may be expected and are very well correlated 
with the sigma-t distribution, L' , M* and N' in Figure 16. 
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correspond to 10 -me ter level. 
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FIGURE 11. Objective analysis of sea surface temperature, California coast, for September 1969 
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FIGURE 12. Objective analysis of temperature at the 10-meter level, California coast, for 
September 1969 
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FIGURE 13, Geographic position of the grid for the Low Latitude Analysis 
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FIGURE 14. Objective analysis of sea surface temperature for low Latitude. 
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FIGURE 15. Objective analysis of sea surface salinity for low Latitude. 
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FIGURE 16. Objective analysis of sea surface sigma-T for low latitude. 



V. 



r.ONr.T.TIST0NS AND SUGGESTION FOR FUTURE RESEARCH 



A. CONCLUSIONS 

A small scale, high resolution and low computing time objective 
analysis model has been developed. More tests and field verification 
are required to show all of its advantages and shortcomings but tests 
described and interpreted in this study suggest that it has great 
possibilities to become a basic step in forecasting oceanographic 
parameters on a tactical unit level. One desirable feature of the 
results was the detail obtained in contrast with that obtained in a 
subjective analysis. The required input and simplicity enables it to 
be used to obtain an almost instantaneous analysis during an oceanographic 
survey and therefore to investigate interesting features not foreseen 
in planning stages. Therefore, it could be applied for the purposes of 
reducing costs and increasing scientific results of an oceanographic 



expedition. 

The procedures described need to be refined in order to allow for 
convergence or divergence whenever these phenomena become important. 
However, by understanding how the analytical model (section B) handles 
these features, meaningful interpretations can be made even for these 
features in the results. The data should be corrected for diurnal 
effects but there is not a satisfactory procedure for this (Technical 
report, ASWEPS No. 13). An alternative solution, observations at the 
same hour of the day, may be undertaken by XBT drops. This is perhaps 
unfeasible at most oceanographic stations due to operational costs. 
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B. SUGGESTIONS FOR FUTURE RESEARCH 

A suggestion for future investigation would be to integrate this 

scheme into a forecast procedure. The following steps could be applied 

using the present form of the scheme: 

a. Objective analysis using horizontal grid at standard dep 

b. Conversion to a three-dimensional grid using the vertical grid 

as a complementary tool. 

c. Repetition of procedures a. and b. at equally spaced time steps. 

d. Determination of Lagrangian time rate of change for identificable. 
structures (assuming that some of them are persistent in relation to the 
time step selected as seems to be the case in the neighborhood of 

constant currents)* 

e. Determination of changes in the heat budget for each character- 
istic structure and (relating) it to easily measured parameters such 
amount of clouds in each standard level, sun’s declination (season) and 

hour of the- day. 

f. Measurement of advection by means of changes in. the structures 

in point d. and (relating) it to permanent currents, wind speed, direction 

and duration. 

g. Forecast for the next time step. 

h. Forecast verification and obtention of corrective cefficients. 

In the delineation of the steps leading to a forecast it has been 

assumed that some features may have a relative long duration or per 
sistence but this situation must be investigated for each application. 

If it is proved, identification of persistent or semipersistent features, 
their changes in time and parameters affecting them may be a simple 
approach for detailed forecasts, even for handling tr. sient thermoc 
now practically impossible to forecast. 



55 



PROGRAM OUTPUT 



MESH 

2 



MINIMUM 
6. COO 



MAXIMUM 
14. GOO 



MEAN LAT 
37.000 



28 April - A May 



ROTATION PARAM GRID 
47.00 1 1 



INT 

10 



INITIAL 

8.000 



INCREMENT 

0.5C0 



EPSILON 

0.001 



POSITION CF THE GRID REFERENCE POINT 
LATITUDE longitude 

37 43.200 124 6.000 

OCEANOGRAPHIC STATIONS 



LATITUDE 

36 
26 
36 
36 
36 
36 
36 
36 
36 
36 
36 
36 
36 
36 
36 
36 
36 
36 
36 
36 
36 
36 
- 36 

36 
36 
36 
36 
36 
36 

36 

37 
37 
37 

36 

37 
37 
37 
36 
36 
36 

36 

37 
37 
37 
37 
37 
37 
37 
37 



LONGITUDE 



40.50 
46. 50 

47 . 50 
4 9 o 6 V 
52. 9U 
54. 1C 
55.7 0 
53o 50 
51 o 80 
56. 20 
54. 80 
47.10 
45o 30 
42c 60 

41.20 
39 o 40 

38.20 
36. 00 

33.60 
3 1 o 3 ' 0 

35.30 
39. 7U 
42.70 

47.00 

50.00 

51. 60 

53.30 
54o 30 
55. 70 
5 5.90 

0.40 
3. 8U 
4.20 
57. 40 
2.70 

I. 50 
0. 1C 

59.40 
56.3 0 
53. 20 
5 8.40 
2.60 
10.00 
10.00 
11.00 
10.90 
11.00 

I I . 00 
10.40 



121 

121 

121 

121 

121 

121 

121 

121 

121 

121 

121 

121 

121 

121 

121 

121 

121 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 



53.40 

54.30 
54. 30 
54.50 
54. 10 

54.40 
54 o 4u 
59. 4C 
59. 10 

59.30 
59. 10 

58.40 

58.70 
59. 40 

59.20 
59. 20 
59. 10 

1.60 
4. 60 

7.00 

10.70 

15.20 
1 8 o 70 
14. 70 
10 . 20 

9.70 
9. 20 

7.50 
6.30 
4. 8U 

14.3V 

17.9V 

19.40 
1 0 o Ov 

2 J« 50 
21. 80 
23 .30 

24.00 

28. ov 

33.40 
40. 5L 

48. 40 
. 4 5. 8V 

41.80 

38.50 

3 5 « 2 J 
33.20 

29. 1C 

26.70 



DEPTH 

0.0 
0.0 
0. v 
0.0 
c.o 
0.0 
0.0 
1.0C 
G.O 
0.0 
c.o 

G.O 
0.0 
0 . 0 
0.0 
c.o 

\J 3 0 
0.0 
o.u 

0 o 0 
0 .0 
U.U 
0 .0 
0.0 
0.0 
c.o 
c.o 
1.00 
o.u 

0 o 0 
Coo' 
0.0 
c.o 
c.o 
0.0 
0 . 0 
0 .V 
1.00 
0 o 0 

G.O 

0.0 
0.0 
G o 

c.o 

0.0 

C o L 4 

c.o 

0.0 

0.0 



TEMP SOUND VELOC 



11.64 
10. 28 
10. 31 
10. 62 
11.31 
1 1.07 
11.02 
10.27 
9.91 
10.82 
10.53 

9. 93 
10.03 
10® 09 
1C. 4 3 
1 0. 60 
10.7 5 

1 i x~i 
11. v4 
1 0 ® 0 ( 

10 .39 

10.41 
9.9b 

10. 12 
10.53 
1C . 25 
1 0 o 65 
10.6 6 

1 1 . 40 
11.24 

9. 62 
9. 5- 
9.49 

11.42 
9.38 
9.56 
9 o 7' 1 
9.85 

10.41 
10.48 
12. C4 
12.45 
12 . 22 
12. 20 
11.70 
11.21 
10.81 
10.21 

9. 85 



1493.20 

1488.50 
1488.70 

1489.80 

1492.10 

1491.20 

1491. 10 

1488.70 
1487.4C 
1490.30 

1489.40 

1487.50 

1487. 80 
1487.90 

1489.10 

1489.70 
1 4 9 C • 5-.j 

1492.40 
1491 . 20 

1487.60 
I486® 3u 

1488.70 
1 4 8 7 . 4 • J 

1487.70 

1489.40 

1488.30 
1489. vG 

1490.60 

1492.70 
1491.90 

1487. 30 

1486.70 
1486. 50 
1 48-feo 30 
1486. 3U 
I486. 9C 
1487.40 
1 4 8 7 . 9 0 

1489.60 

1489.70 
1494. 4C 

1495.70 
i 4 9 5 o 0 V 
1495o 3v 
1493o3i.- 

1492. 60 

1489. 90 
1489.20 

1487.90 



G.O 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
G .0 
0.0 
6 c 50 

12.90 

2 ( o <0 

29.70 

22.30 

2C.20 

24.20 

"i AG 

' 5o 62 

12.60 

64.90 

53.90 
48.50 

6 3 • G 0 

7 1.70 
79.10 

54.00 

2 6 o 5 0 

9.70 
5.10 
5.7 0 

21.40 
22 » 5o 
23. UO 
7.60 

36.00 

3 3.20 
26.50 
19.30 

v ® 0 
18. 4V 

65.90 

62.90 
t>6 . 50 
62. 5C 

62.70 
9. 50 
9 o 6u 

10. 2C 
16.20 
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37 

37 

37 

37 

37 

37 

37 

37 

37 

37 

37 

37 

37 

37 

37 

37 

37 

37 

37 

37 

37 

37 

37 

37 

37 

37 



10.90 


122 


2 5.20 


1.00 


9.25 


1486.00 


26.80 


15 . 80 


122 


27 . 50 


C • 0 


10.26 


1489.40 


11.00 


21.60 


122 


31.70 


CoC 


10 . 54 


1 490. 20 


9 . 50 


31 . 90 


122 


33 . 4 l > 


c • o 


9.92 


1494.30 


19.80 


31.60 


122 


35 . 20 


c.o 


11.12 


1492. 30 


4.10 


32 . 10 


122 


37 . 10 


c.o 


11.65 


1493. 60 


5.00 


31.90 


122 


3 8 . 70 


c.o 


11.35 


1492 ® 3 u 


5.30 


32.20 


122 


41.30 


c.o 


11.41 


1492.60 


22.10 


31.90 


122 


43 . 10 


c.o 


11.18 


1492.50 


31.80 


31.60 


122 


46.90 


c.o 


9.76 


1486.20 


62.30 


32.00 


122 


49 . 90 


0.0 


10 . 49 


1469.60 


63.10 


32.30 


122 


54.50 


c.o 


10.98 


1491.20 


61.20 


31 . 2 C 


122 


58.00 


0.0 


lx . 88 


1493.90 


68.20 


32.30 


123 


3 . 80 


0.0 


12 . 27 


1495. 10 


72.60 


32 • 40 


123 


9 . 30 


c.o 


12.06 


1494.40 


63.80 


36.80 


123 


15.40 


c.o 


11.07 


1491.60 


64.70 


42.20 


123 


7.70 


c.o 


10.46 


1489.40 


77.20 


47 . 70 


122 


59 . 70 


0.0 


10.19 


1489.20 


6.10 


52.60 


122 


52 . 5 C 


c.o 


9.94 


I 486. 30 


12.00 


53 . 2 C 


122 


49.20 


0.0 


9.94 


1486.30 


9.80 


49.00 


122 


48 . 8 C 


C.o 


9 . 58 


1467.00 


17.70 


46.40 


122 


46.90 


0.0 


8.95 


1484.60 


49.60 


43 . 50 


122 


44 . 50 


CoC 


10.71 


1490.30 


16 . 0 C 


4 -j . 20 


122 


42 . 10 


0.0 


10.81 


1490.00 


19, 10 


37 . 50 


122 


40.20 


0.0 


10.80 


1490.30 


23.00 


37 . 30 


123 


36 . 60 


c.o 


10 . 62 


1469.50 


18.10 
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relative 



118. C91 
113.213 
112.482 
110.691 
108. £42 
1C7.601 
106.421 
105.316 
106.722 
103.396 
104.529 
110.541 
111.693 
113.286 
114.419 
115.735 
116.667 
116.914 
117.025 
117.409 
112.468 
106.799 
102. 698 
101.723 
1C1.991 
101. C94 
100.123 
100.318 
99.948 
100.619 
92.153 
87.7C5 
86. 595 
96.689 
87.093 
87.262 
87.469 
87.599 
87.360 
87.012 
79.341 
71.966 
67.971 
70.150 
71.217 
73.088 
74.104 
76.338 
78.084 
78.536 
73.7CC 
67.170 
58.712 
57.951 
56.550 
55.825 
54.189 

53.428 
51.577 
49.650 
46.924 
45. 822 
41.858 
38.788 



POSITION IN THE GRID 

y PARAMETER 



34.677 
38.244 
38.926 
40.378 
42.727 
43.370 
44.461 
40.040 
39.C56 
41.940 
41.102 
26.259 
34.856 
32. 605 
31.767 
30.539 
29,779 
26.818 

23.429 
20 c 4 59 
21.C26 
21.399 
21.401 
26.671 
31.345 
32,729 
.s4 * 18u 
35.855 
37.511 
38.523 
3 6 5 0 A A 
36o 261 
35. 653 
36o 51(3 
33.992 
32.414 
30 . 5 
29.697 
24o 896 
19.978 
19,378 
17.629 
24.195 
26.531 
29.141 
31 . COG 
32.236 
34. 631 
35. 623 
36.340 
38. 839 
40.343 
40 o 376 
45. 120 
44.351 
43.280 
41.966 
40.710 
38. 286 
36. 807 
34.325 

31. 530 
28.893 
25.749 



11.640 
1C. 280 
1C. 3 10 
10. 620 
11 . 3 1C 
11.070 
11.02C 
10.270 
9. 9 1C 
10.820 

10.530 
9.930 
1C.03C 
10.090 

10.430 
10. 600 
10.750 
11.370 
11.040 
10. 070 
10.390 
10.410 
9.980 
10. 120 

10. 530 
10. 250 
10. 600 
10. 860 
11.400 
11.240 
9.620 
9.550 

9.490 
11.420 

9.3 80 
9. 560 
9.710 
9.S5C 
1C. 4 It 
10.480 
I2.u40 
12.450 
12. 220 
12.300 
11.700 
11. 210 
10.810 
1C. 210 
9.8 50 
9.3 50 
10.260 
10. 549 
9.920 
11.120 
11.650 
11.350 
11. 410 
11.180 
9.769 

10. 490 
10.980 
11.880 
12.270 
12.060 
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32.247 
32.494 
3 2 o 820 
33.170 
34.529 
37.818 
40.754 
44.183 
47.903 
50.913 
20.322 



25. 188 
33.368 
41.792 
49.340 
51.676 
49.C45 
48*33.1 
47.605 
46.956 
46 . 224 
13.147 



11.070 
10.460 
10. 190 
9. 940 
9.940 
9. 5 80 
8.950 
lu.710 
10. 810 
10.800 
10.620 



I 
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INTEPPOLATED 



8*00000 
16*00000 
24*00000 

32.00000 
40. C 0000 

48.00000 
56.CC0CC 
64.00300 

72.00000 
80.CC0C0 
88* OCCoO 
96*00000 

104.00O00 
112*00000 
120.0Q0u0 
128*00003 
8.CG000 
16*00000 

24.00003 

32.00000 
40. coy oo 
48* 0 0 0 w u 
56. CO 300 
64 « 0000 0 

72.00000 

80.00003 

68.00000 

96.00000 

104.00000 
112.C0CG0 
12C.CC CUp 
1 2 8 . C 0 0 0 3 

8. C 0000 

16.000 0 • j 
24. CO 000 

32.00000 

40.00000 
48.00000' 
56.GCCCC 
64.00030 
72. C CLOU 
80. OCC00 . 
88.0CG0O 
96. C COCO 

104.0C0O-C 
1 12.0 COL c 
120.0L OvO 
1 2 8. L COCO 
8.CCC00 

16.00000 
24.LG0C0 
32.00030 

40.00000 
48*00003 
56 . C COCO 

64.00000 
72.CC00C 
80. C COCO 
88.0C0CC 
56. COO 30 

1C4.000GC 
1 12. CGOCG 
120.00000 
128.00CC0 



VALUES FCR THE GRID 
Y PARAMETER 



8. COO 00 

8.00000 
8 . UOOGO 

8.00000 
8. CCOGO 

8.00000 
8. COG 90 
8.0C000 
8.C000C 
8.CG0C0 
8 o CO 000 
8.00CCC 
8. CCL 00 

8 .00000 
8 o Ci OOo 

8. 1 LOCO 
16.COOOC 
16.GC 0 CO 
16.00300 

16.00000 
1 6 ® 0 C C 0 0 

16.00000 
16o COL 03 
16.0 CL 00 
16 « C 00 00 
1 6 . C v 0 00 
16. COO 00 
16. CCD CO 
16. COC 30 
16.C0CCC 

1 A \ . * >Q ' i Ly 
t 7 V 7 , . ,< 

lOol ‘ J O v3 V 

24o COOOO 
24. COO 00 
24.OCC0 3 
24 o COC C 3 
24. t COCO 
24.0000'.- 
24. CvCOO 
24.UC000 
24.03000 
24. LLOOC 

2 4 o 00 0 0 0 

24. C 0000 
2 4 o v v. 0 0 0 
24. 00 GOO 
24.CC00C 
24.CCGCU 
32 . CjG 00 
52. C LOCO 
32. LOG 00 
32. COL GO 
32. C 3000 
32 o UuC 03 
32. COO CO 
32.CC DOC 
32 « COC DO 
32 • COL DC 
32.0300*3 
32. COC DC 

32. LGOjG 
32. GOO CO 
32.CC COO 
32.000 CO 



0.0 

0.0 

0.0 

0.0 

0.0 

o.c 

0.0 

o.c 

o.c 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

1 0. 620uu 
1C. 6 2 u v> 0 
0.0 
0.0 
0.0 
c.o 
0.0 

12. 453 3 0 
12. 04000 
10. 48300 
0.0 

0.0^ „ 
8*95623 

0. e 

n . f' 

V W v 

0.0 

o. c 
o.o 

1 1 o 07CG0 
12.06000 
0.0 
0 o 0 

12.220CU 

12.3DL00 

O.C 

10.41C0C 

0.0 

10.^7425 
10. 39000 
10. 9 1 o 5 2 
U.O 
0.0 
0. 0 

c.o c , 

2.05547 
12. 27000 
1 1.43289 
0.0 

11.14965 
11.32377 
9.05577 
9.63905 
9.92117 
10.33595 
9.84487 
12. 34712 
G * 0 
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8.0CCCC 
16* COOGC 
24. CCOOO 

32.00000 
4 0 • 0 \j 0 0 0 

48.00000 
56.0C0C0 

64.00000 
72.CCCC0 
80.0C0C0 

88.00000 

96. OCUOO 

lC4«COOvO 

1 12.00000 
120 . G0CC0 
128.00030 

8. OttOCO 

16.00000 
24. CCOOO 
32.CCOOC 
40.000C0 

48.00000 
66.CC00U 
64 • COOuO 
72. CCOOO 
80. C OoOO 
88.00000 
96.CCCC0 
104.COOOO 
112.CCCOO 
1 2 0 . 0 v 0 v • C 
128.CCCC0 



40. 
40. 
40. 
40. 
40 « 
40. 
40. 
40. 
43. 
40. 
AJ e 
4G . 
40. 
40. 
40. 
40 . 
48. 
48. 
48. 
48. 
48. 
48. 
48. 
48. 
48. 
48. 
48. 
48. 
48. 
48. 
48. 
48. 



CCOOO 

CCOOO 



CCOOO 

ccooc 

CCOOO 
00000 
000 oc 
000 00 
CCOOO 
COO 00 
COG 00 
CCOOO 
00000 
coe go 

COCOO 

cocuo 

CGCCU 

CCOOO 

00000 

CCOOO 

CuGOO 

CCOOO 

CCOOO 

C 0 0 0 u 
00000 
COO 00 
oeooo 
COOuG 



V s Kj 'J 



_ - . 00 
CCOOO 

v»uOoO 

CCOOO 



0.0 

0.0 

0.0 

10. 19000 

0.0 

10.26933 
11.43151 
10. 5400C 
10, 26000 
9. 3 5000 
9. 55000 
0.0 

10.54430 
10.45130 
0.0 
0.0 
0.0 
0 o 0 
0.0 

10.00981 
8. 11982 
10.72071 
8. 92048 
10. 39016 
0.0 
0 . 0 
0.0 
0 « 0 

11.02300 
12. 20851 
0.0 
C.O 
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U/K 


V/K 


RATIO 


K 


L 


0* 6C0 


0.0 


0.600 


10 


2 


G.6GC 


O.G 


C. 600 


11 


2 


0.CC6 


-0.600 


0. 600 


12 


2 


Co 6C0 


0.0 


0. 600 


13 


2 


0.600 


OoC 


Oo 60 0 


14 


2 


0.6CC 


OoC 


0.600 


9 


3 


0.6CC 


0.0 


0.600 


10 


3 


0.591 


-0.102 


0.600 


11 


3 


C.6C0 


OoO 


0.600 


12 


3 


-0.053 


0. 599 


0. 600 


13 


3 


0.600 


O.C 


0.600 


14 


3 


0.6C0 


0.0 


0. 60C 


15 


3 


C.6CC 


0.0 


0. 600 


7 


4 


-0.374 


0.469 


0. 600 


8 


4 


0. 6C0 


O.C 


0,600 


9 


4 


0.6CO 


0.0 


0. 600 


10 


4 


-0.155 


Co 580 


Co 60Q 


11 


4 


0.152 


-0. 569 


0.600 


12 


4 


0.600 


0.0 


Co 600 


13 


4 


0.6(0 


O.C 


0.600 


14 


4 


0, 6 CO 


c.o 


0.600 


5 


5 


0.610 


C e C 


0. 600 


6 


5 


C.6CC 


OoC 


0 o 600 


7 


5 


0.547 


— 0 o 246 


C. 600 


6 


5 


0.546 


-Co 249 


0 . 6 u 0 


9 


5 


0.600 


O.C 


C. 600 


10 


5 


0.6(0 


OoO 


0 o 6 v 0 


11 


5 


0.600 


0.0 


0. 600 


12 


5 


-0.064 


Oo 697 


0.6 GO 


13 


5 


0.6CC 


0.0 


0 o 6 G v 


14 


5 


0.600 


0.0 


0* oCO 


15 


5 


0. 6GG 


O.C 


0. 60 o 


8 


6 


0.600 


OoC 


0. 600 


9 


6 


C . 600 


O.C 


0 e 6>> 0 


1C 


6 


0.6 CO 


0 o Vi 


0.600 


i i 


6 


0.600 


O.C 


Oo 6o0 


12 


6 


C. 6C0 


OoC 


0.600 


13 


6 


0.6CC 


0.0 


0. 600 


14 


6 
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FINAL VALUES FOR GRID POINTS 



TEMP 

9.69 
9 . 89 

9.69 

9.89 

9.89 

9.89 

9.69 
9 . 89 

9.89 
9 . 64 
8.79 

11.05 

12.88 

11.46 

9.89 
9 . 89 

9.89 

10 o 62 
10.62 
9 . 89 

9.89 

9.89 

9.89 
9 . 89 
12.45 
12.04 
10.48 
10.66 
10.81 
8.96 
7 . 98 

9.89 

9.89 

9.69 
9 . 89 
11.07 

12.06 

9.89 

1 C . 61 
12.22 
12.30 

10.89 
1 C . 4 1 

10.34 

10.47 
10.39 
10.92 

12.90 

9.89 

9.89 
9 . 89 

9.89 

12.27 

11.43 

10.75 

11.15 

11.02 

9.06 

9.64 

9.92 

10.34 
9 . 84 

12.35 

9.89 



L 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 



K 

1 

2 

3 

4 

5 

6 

7 

6 

9 

10 

1 1 

12 

13 

14 

15 

16 

1 

2 

3 

4 

5 

6 

7 

6 

9 

10 

11 

12 

13 

14 

*i s; 

16 

1 

2 

3 

4 

£ 

6 

7 

8 

9 

10 

1 1 

12 

1 3 

14 

1 5 

16 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

1 1 

12 

13 

14 

1 5 

16 
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D 




9.89 
9. 89 

9.89 
11'. 19 
7.51 
10.27 
11*43 
10.54 
10.26 
9. 35 
9.55 
1C. 38 
10.54 
10.45 
9.30 
6.45 
9. 89 

9.89 

9.89 
1C. 01 

8 o i 2 

10.72 
8.92 
10.39 
1C. 10 
8.57 
11.01 
10.13 
11.02 
12.21 
6.03 

9.89 



5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 



1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

1 1 

12 

13 

1 4 

15 

16 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 
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program output 



MESH 

2 



MINIMUM 

6.000 



MAXIMUM 

14.000 



MEAN LAT 
37.000 



INT 

10 



INITIAL 
80 COO 



28 April - 4 May 



ROTATION PARAM GRID 
47 oOO 1 1 

POSITION OF THE GRID REFERENCE POINT 

LATITUDE LONGITUDE 

37 43o 200 124 60 COO 

OCEANOGRAPHIC STATIONS 



INCREMENT 

0.500 



EPS ILON 
0.C01 



LATITUDE 

36 

36 

36 

36 

36 

36 

36 

36 

36 

36 

36 

36 

36 

?6 

36 

36 

36 

36 

36 

36 

36 

36 

36 

36 

36 

36 

36 

36 

36 

36 

36 

36 

36 

36 

36 

37 
37 
37 
37 
37 
37 
37 
36 
36 
36 

36 

37 
37 
37 



40.50 
42. 20 

43.80 
45. ON 

46, 50 

47. 50 
4 9 o 8 L 

51.40 
52.90 
54. 10 

55.70 
56 . 20 
54. 80 
53. 50 

51. 80 
47. 10 
45 o 3^ 
44. 30 
42 o 60 
41. 20 

39.40 

38. -20 

36.00 
33. 6 G 
31. 30 
3 5.30 

39. VC 

42.70 

47.00 

50.00 

51.60 
53. 30 

54.30 
5 5. 7 0 
55.90 

0.4 r 

3. 80 

4. 20 
2.70 

31.40 
33.1C 

0 . 10 

59. 40 

56.30 
53. 20 

58.40 

32.60 
10. CO 

10.00 



LONGITUDE 

121 
121 
121 
121 
121 
121 
121 
121 
121 
121 
121 
121 
121 
121 
121 
121 
121 
121 
121 
121 
121 
121 
122 
122 
122 
122 
122 
122 
122 
122 
122 
122 
122 
122 
122 
122 
122 
122 
122 
122 
122 
122 
122 
122 
122 
122 
122 
122 
122 



10.70 

15.20 
18. 70 

14.70 

10.20 
9. 7C 
9. 20 
7. 50 
6.30 

4 o 80 

14.30 
17. 90 

19.40 
20.50 
21 . 80 

23.30 
23. 3C 
24. 00 
28. 60 
33. 4^ 
40. 50 

48.40 

45. 80 
41.80 



DEPTH TEMP SOUND VELOC 

1490.00 
1469.40 

1487.00 

1 483.80 
1488.60 

1487.80 

1489.80 

1491.30 

1491.90 

1491. 50 

1491.20 

1489. 90 

1469.20 

1487. 30 

1486. 90 
1 486 ,4*j 
1 4 6 . 3 w 

1467. 20 
1489.10 
1499. 2 r ' 

1489.90 
1488. 6C 

1488.80 

1491.40 
1 4 8 7 o 60 

1488. 50 

1488.90 

1487.40 

1486.90 
1485o 8‘J 

1486.70 

1489.20 

1465.40 
1487. 10 

1488. 70 

1466. 50 

1486.70 

1486 . 20 

1486.40 
1486.00 
1487. 10 
1487. 10 

1488.20 
1489.60 

1489.80 

1493. 20 

1495.80 

1475.20 

1494. 70 



53. 7C 


10.00 


1C. 61 


54. 7C 


11.00 


10.51 


54. 60 


11.00 


9.76 


54. 30 


l n .0C 


1C. 28 


54 . 30 


9 0 0 0 


10.27 


54. 30 


11.00 


10.03 


54 0 5u 


9.00 


10.58 


54. 30 


10.00 


11.03 


5^.10 


11.00 


1 1 0 1 8 


54.40 


9.00 


11.08 


54. 40 


9.0C 


1 1 0 0 1 


5°, 30 


9.00 


10.63 


59. 10 


lr.oCO 


10.42 


59. 40 


10 0 G 0 


9 0 85 


59. 10 


10,00 


9.78 


58. 40 


10.00 


9o 6 :, 


58.70 


10.00 


9 » 6 r 


69. 10 


11.00 


9 0 8 7 


59. 40 


10. oc 


10.08 


59. 2C 


11.00 


10.41 


59.20 


12.00 


1 n . 63 


59, 10 


10. CO 


10o 51 


1. 60 


11.00 


lOo 47 


4. 60 


8.00 


11.07 


7.00 


10.00 


10. or 



8.00 

10.00 

10.00 

10.00 

7.00 

8o^0 

9.0C 
11.00 
11.00 
1C. 00 
10.00 

8.00 

7.00 

8.00 

9.00 
9.00 
9 of. n 

11.00 

8.00 
9.00 
9.00 

10.00 
10.00 
9. CO 



10.45 
9,98 
9.89 
9.48 
9.77 

10.14 
9. 40 
9.70 

10. 3c 
9.32 
9.47 
9. 29 
9.36 
9 ® 2 Q 
9.60 
9 o 60 
9. 84 

10.28 

10.46 

11. 57 

12,45 

12.23 

12. 1 ) 



0.0 

0.0 

OoG 

c.o 

A A 

r O J 

0.0 

0.0 

C.O 

C.O 

c.o 

0.0 

^.0 

0.0 

c.o 

0.0 

a 

L ft V. 

*# © w 

0 o 0 

c.o 
0.0 
c.o 
0.0 
c.o 
0 .0 
0.0 

0.0 

C n 
© -> 

0.0 

2.0 

0.0 

c.o 

0 o 0 

c.o 

0 o 0 
0 o 0 
0.0 
0.0 

Q.° 

c.o 

0.0 

OoC 

0.0 

0.0 

0.0 

0.0 

0.0 



0.0 

0.0 
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122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

123 

123 

123 

123 

122 

122 

122 

122 

122 

122 

122 



38o5C 
35o 20 
33. 20 

31.30 
29o 10 

26.70 

25.20 
27. 50 
29. 60 
3 1 o 70 
33. AC 

35.20 
37. 1C 

38.70 

41. 30 

43. 10 

46.90 

49.90 
54. 50 
58.00 

3.30 

9.30 
15.40 

7.70 

59.70 
52. 50 
49. 20 
48. 80 
46. 90 

44. 50 
42. 10 



9.00 
10.00 
11.00 
10.00 
9.00 
10.00 
10. oc 
9. CO 



8 o CO* 

9.00 

3.00 
11.00 
11.00 

°o 00 

10. oc 

9.00 
9.00 
9, CO 
9.00 
8 o OC 

10.00 
11.00 
10. oc 
10.00 
8. CO 
10.00 
9.00 
9 .00 
10.00 



11.18 
10.27 
10.57 
9. 1C 
1C. 09 
9. 49 
9. 32 



10 

9. 

9. 

8 . 

9. 

9. 

9. 

9. 



05 

87 

51 
98 
32 
78 

52 



10. 59 
9.75 
10.46 
10.91 
11.04 
12. 25 
11.98 
10.61 



9.36 

9.65 

9.95 

9.92 



9.44 
8.93 
10.22 
9. 82 



1491.80 
1489. 30 
1490. 50 

1465.80 

1488.90 

1436.60 
1486.20 

1488.80 
1488. 1C 

1486.70 

1485. 60 

1465.90 
1487.10 

1486.80 
1487.40 

1490.60 

1487.40 

1489. 70 

1491.30 

1491.40 

1495.30 

1494.30 

1490.30 
1485o 40 

1487.30 

1488.40 

1488.40 

1 436.40 

1484.60 
1 4 38 . 8C 

1487.90 



v • O 

c.o 
0.0 
c.o 
0.0 
0.0 
c.o 
0.0 
C o G 

o.c 

0 o 0 

O.o 

o.o 

c.o 

c.o 

o.c 

0.0 

o.c 

c.o 

0.0 

c.o 



0.0 

0.0 

OoO 

o.c 

c.o 

0.0 



o n 

’J o - 

0 o 0 
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FINAL VALUES FOR GRID POINTS 



TFMP 
9 o 70 
9 o 70 
9 0 70 

9.70 
9 o 70 
9o 70 
9 0 70 

9.70 
13 o 42 
lie 84 
ICo 32 

8.79 
7 o 62 
8o09 

8 o 28 

9o 70 
9 0 70 
9 0 7 0 
9o 70 

9.70 

9.70 
12.77 

9.70 

9.70 
12.83 
11. 57 
10.46 
9.65 
9. 32 
9. 23 
11.75 
9 „ 70 

9.70 
9. 70 
6 . 10 
10.61 
11.98 
1C.78 
9.03 
12.23 
12.10 
11.(9 
10.28 

9.98 
10.40 
10.38 
10.55 

9. 51 

9.70 

9.70 

9.70 
9. 36 
12.25 
9. 57 

9.80 
10. 29 
10.05 
10.42 

9.59 
9. 59 
9. 81 
9.61 
10.54 

6.98 



L 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

5 

2 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 



K 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

1 1 

12 

13 

14 

15 

16 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

1 1 

12 

13 

14 

15 

16 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

1 5 

16 
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9.70 

9.70 
11.61 
9.65 
7.67 
9.38 
10.20 
9.86 
9.92 
9.32 
9. 47 
10.14 

9.89 
9.86 
7.20 
8.26 

9.70 

9.70 

9.70 

9.90 
8.26 

10.19 
9.76 
9.14 
10.46 
7.49 
10.63 
7.64 
11.01 
12 . 64 
8 . 10 
9 . 70 



5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 



1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

1 6 
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PROGRAM FOP OBJECTIVE ANALYSIS OF CONSERVATIVE PARAMETERS 
OF THE OCEAN 



REAL*8 TJTLE 
LTG 



COMMON/OBS/ XOI 120) t Y0( 120) t 0BS1 1 12C I 

CCMMON/OATA/LAtV lor') 'iOMINL }00 I tLONGI 100 I , OMING (100 1 1 

DERJHI 1001 ig?i < IK > iVYk > 



2 DEPTH ( IOC ) , UB1 l 1UC j ,un<c i ^\’’X Y UJl 

COMMON/PARAM/X(20) , V 2J) 

DIMENSION CL( 10 ) , LTG(3 ) ,TI L - ( 1 2 ) 



C 

C 

C 

C 

C 

c 

c 



PROGSlS^Sl^I^FpJ 5 ANALYSIS ^ d § n F or A missi ng L data?var?able 
GRIO^SIZE°AN ROTATION of GR I Do X SC THE ATI CNS C AVA I L ABL 

ROTATION OF THE GRID* 

READ (5,898) 

READ (5,329) 

RICK IS THF MEAN 

WRItI{!i326) MESH.P.MIN.RMAX 

8UBI?:S1ti ANG^lf AS, MI, INT, SINT, DELTA, EPSI 

wr t t P ( 6 ! ^3>T Sang, i p a r , i m » i nt , sint , delta, etsi 

READ ( 5 , 2®2 ) IL AT , RM INL » 1 LONG , RM I NG 
WRITE (6, 2 9-3) 

WRITE ( 6 » 293 ) I L AT , RMI NL , I LONG , RM ING 

READ ( 5 , 280 ) XSC , YSC , NMR , MMR , NR 

THET=0o01745*ANG 

RICK=0o01795*RICK 

JSC = 1 

KJ1^9 

READ (5, 272) ( LAT ( I ) , OM INL ( I ) i, LONG ( I ) , OM I NG ( I ) , DEPTHt I 

2 , OBI ( I ) , 0B2 ( I ) , 0B3 ( I ) , 1-1 »NR ) 

WR ITE ( 6 , 270 ) 

WRITE (it 279) ( LA T( I ) .OMINUI ) , LONG ( I ) , OMI NG( I ) , DE PTH( I 



THE 



(TITLED ,1=1,12) 

MESH, RMIN, RMAX, RICK 
LATITUDE FOR THE GRID 

RICK 



C 

C 

C 

c 

c 



TO OblAiN \ r tr ncrucuis 

SIS^fcHICAL^OIJTION IS 

IE DATA CARD 
OBSERVATIONS FOR EACH S T AT I ON ON 
HORIZONTAL GRID MI MUST BE ONEo 
WRITE (6,321) 

WRIT E ( 6 , 322 ) 

I F ( t LAT-LAT (^ J ) ) 404,406,*-05 
LAT(J)=LAT(J)-1 ^ ^ 

CMINL(J)=0MINL(J)+6uo0 

GO TO 903 , , 

LAT ( J ) = L AT ( J ) + l „ . 

CMI NL ( J ) =OMI NL ( J ) -60o 0 

?F (1 LONG-LONG ( J ) )9D7,909,908 
LGNG( J)=L0NG(J)+1 
CMINGt J ) =OMI NG ( J )-60, 



903 

909 



9C5 



906 

908 



,0 



9 C 6 



9C7 



GO TO . , 

LONG( J ) = LONG ( J ) - 1 A _ 
CMI NG ( J ) = 0 M I NG ( J ) +6O0C 

GO TO 906 
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409 X0 ( J ) =( RMING-OMI MG ( J) ) *COS ( RI CK ) 

XOREL=XO ( jV=COS^THET )-Y0 (J)*SIN(THET) 
YOREL = XO(J) ; ' : SIN(THET)+YO(J)* ; COS(THET) 

IF(MIcEQ.l) GO TO 411 

TKE *i8 5 2.4 

Y0( J)=DEPTH( J) 

GG TO ^l r ' 

411 X0( J ) =XCREL 

YO ( J ) =Y OREL #i/ /w 

410 IF(IPAR-2)412,414,416 

412 OBS 1 ( J ) =0B1 ( J ) 

GO TO 418 

414 OBS 1 ( J ) =0B2 ( J ) 

GO TO 418 

418 WRITE (*6, 32 3) X0( J ) , Y0( J) , OBS 1 ( J ) 

419 CONTINUE 

2 XC=XSC / 2o 0 

YC=Y SC / 2o 0 c 

XR=(XC**2+YC**2 ) **0o 5 
WRITE ( 6 , 349) 

WRITE ( 6 , 322 ) 

DO 15 L=1 , MMR 
YG=YSC*L 

4 DO 15 K = 1 » NMR 
XG=XSC*K 

N= 1 

DO 8 1=1, NR 
XOB=ABS ( XO ( I )-XG) 

IF( XOBoGT. XC ) GO TO 8 
YOB=ABS ( YO ( I )-YG) 

I F { YOBoGToYC) GO TO 8 
7 (N)--OBSI ( I '< 

X( N) =Xu ( i J — XG 
Y( N) =Y0( I ) -YG 

5 N=N + 1 

8 CONTINUE 

IF (N-2) 10,12,14 
10 GVALl(K,L)=CoO 

12 AP=Tx( 1 )**2+Y( i)**2)**0 o 5 
I F ( AP • GT o 7« 0 ) GO TO 10 

IF^GVALUK,L)o LToRMAX) GO TO 15 

GVALl(K,L)=CoO 

GO TO 1 5 

14 N=M-1 

CALL INTERPtN, XR, POLY) 
^feWI^^T.RMAX, GO TO 15 

15 WRITE ( 6 \ 350*) XG, YG, GVALl (K,L) 

N1=NMR-1 

N2=NMR-2 

Ml=fv'MR-l 

R2=MMR-2 

c p NI IS I THE V RATIO F CF THE TWO COMPONENTS OF THE 
RAT 10=0 o 6 
WRITE (6,319) 

DO 150 M=1 1 M J1 
DO 150 L=1 , MMR 
P=0o 27 

RAY=C* ( RAT IO**2~RAX**2 )**0o 5 

41 ??(^ C alK;lTgt 0 o 0 i» GO TO 14" 

I F ( G V AL 1 { K ! L + 1 ) I L T o 0 o C 1 ) GO TO 1104 
IF (GVALl ( K, L+2 ) o LToOoOl ). GO TO 11U4 



VELOCITY VECTO 
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GVAL1 < K, l) = PAY*YSC*(GVALl(K,L + l ) -GV AL 1 ( K , L +2 ) )+2»0- 

GO TO 1102 

I F ( ABS ( GVAL1 (K » L ) ) o GT « RMAX ) GO TO -10 2 
GO TO 149 ^ A 

1102 GVALl(K,L)=0o0 nnrn Tn , 4Q 

1104 I C (GVAL1 (K,L-1 )oLToOovU GO 10 149 

I F ( GVAL 1 ( K » L-2 ) « LTo Co 01 ) GO TO llo6 

GVALl < K,L)=RAY*YSC* (GVALl (K t L-2 ) -GVALl (K,L-1 ) > +2 *°* 

2 lRABMGVALlTKyU } •§{3 aX?G0 TO 1105 

I F { ABS ( GVAL 1 ( K » L ) ) o GT o RMAX ) GO TU 

GO TO 149 

1105 GVALl(KtL)-lo l ^ -yp 149 

ns6 jsiswicHU wiOTK',c- Htf?jw K,Lti ”' 2 * 0-MV * vs c* 

«si 

GO TO 149 
GVALl i K t L ) =, 'o0 
GO TO 149 

I 7 i N M R- ( Kf 1 ) ) r> 1 » 5 4 , 4 ? 

Ir'.LocQwOMR) GO TO 185 



11C8 

50 

51 
54 
42 
80 
83 
81 



GJ T< 
1 r . L . 
If ( L, 
If V L, 
if (L 



Li') u OOR) GO TO 01 
£0,1/ GO TO 11/.. 
GT o HI ) GO i 0 I 
EQ.l) GO TO ilO 



I F ( G V A L 1 U+ 1 ^ L + 1 ) . LT o 0 , 1 ) GO T 0 8 1 
IF(GVALl(K+2»L)oLToGol) GO T 0 8 £ 

I F ( GVAL 1 ( K+ 1 ? L - 1 ) o L I oVj o 1 ) GO Tu 81 



N = 7 



TO 103 



1C 08 
ICO 



103 

104 



I F ( GVAL 1 ( K+l »L lolTofo i ) GO 
I F ( No LT , 7 ) GO TO ICC 

CALL MOLE ( MM»K » L *XSC * Y S C » D IV ) 

GVAtllKiL)=GVALl(K»Ll -RAYM GVAL 1(K+1»L+1 )“GVAL1(K + 1» 

2 1FUBS(GVAL1(K , .L ) j •SToRMAXI GO TO TOGS 
IF(ASS(GVAL1 (K»L) ) ® LT.RMIN) GO TO 10<-S 

GO TO 149 



GO 

GO 



TO 

TO 



103 
10 3 



GO TO 11' 



1C 39 
108 



110 

111 



GVALl ( K i L ) =C .0 
I f( GVAL1 ( K,L~1 ) e LTo % 1 ) 

I F ( LoEQoMMR ) GO TO 104 
If (GVAL 1 (KtL+1 ) oLToOol ) 

N=7 

GO TO 104 

IRKoLT.2) GO TO 132 
I f ( GVAL1 ( K- 1 » L ) , LI .0.1) 

IF (N.LT.7) GO TO 108 

CALL MOLE ( MR ,K , L ♦ XSC ♦ YSC » D I V ) AX 

DI V= ( GV A Li ( K+ 1 , L ) -GV AL 1 ( K-l » L ) , I ' "R A * 

2 K > L+l 1 -GVAL I ( K , L- 1 ) ) 5% ^A / / ( 2o0^ VSC ) 
GVALl ( < t L t -OVA LI (K, U-'J IV 
I - 1 ABS ( GVALl ( K , i. ) ) ob o 9 MAX) 

I F ■. ABS ( GVAL I ( K , L ) ) o L • o Rril N ) 

GO TO 149 



RAX/(2oO~XSC)+(GVALl( 



GO 

GO 



TO 

TO 



1009 

1009 



G JAL 1 ( K j 






av;* 


TP 






1 - ' Lo L V c 




) f 


GO 


A. 


] !' 


>F f GVALl 


t 


K- 


It L 


-1 


1 O 


LT 


N-t 














Go ro -1 














N-0 

l r - 4 o L r o 


± 


i 


f>0 


-f] 


X 


3 2 


i v GVAL A 


t 

*. 


K- 


■ 2 1 L 


) o 


LV 


o ' • 



,0o 1 ) GO TO no 



,1 ) GO TO 114 



75 



I F ' Lo 5 < M* ^ ’ \ ‘0 Ll^ .... ^ 

112 ir\GV«Ll v u- it Lvi )oL u Jo 15 GO TO il 
IF(N.l.T.7) GO TO 113 

CALL MOLE ( MM ? K ? L ? XSC ? YSC ? D I V ) 

GVALl (KtL)^GVALl(K>L) -RAY*? GVAL1 (K-ltL + 1 ) -GV AL1 ( K-l » 
2 IR A8S( GVALljK^L ) ) oGT.RMAX) GO TO 1004 

I F ( AB S ( GVAL 1( K » L ) ) o LT« RMI N) GO TO lLL^t 

GO TO 149 

1004 GVAL1 ( K? L) =0.0 
GO TO 113 

114 N=0 

GO TO 115 
n? N=7 

115 IF(LoEQ.Ml) GO TO 75 

I F ( No LT o 7 ) GO TO 132 rn TO 132 

I F (GVAL 1 ( K -2 ? L +2 ) • L i oOo 1 ) GO TO HI 

N=N+8 

132 IF(GVALHK? L+2 ) I LT 1 1 ) GO TO 134 
IF ( N o LT« 1 5 ) GO TO 136 

CALL MOLE( MM ? K ? L ?XSC?YSC?DIV) 
GfUBMgvlLHK^n.GT.JMAXI GO TO 1003 



1 nADOiyVHLlUNTL f f 

IF(ABS(GVAL1(K,L)) 

GO TO 149 

1003 GVAL1 ( K , L)=0o0 
1 34 IF(K.EQ.NMR) GO TO 75 
GO TO 149 

i36 I F ( Lo GT * Mi } GO TO 7 5 

IPlgvtllWM.foo.l.GO TO 75 
I F ( L« GT • M2 ) GO TO 146 
I F ( K o GT « N2 ) GO TO 148 

I F ( GV AL 1 ( K + 2 » L ) o LT o' : o 1 ) GO TO 148 
I F (GVAL 1 ( K+2 > L +2 ) « LToO o 1 ) GO TO 148 

CALL M0LE(MM,K,L jXSCtYSCtDiy) 

iWtl?K:tV=§vitliK:!: , . _ 2ftIi?52fXv*LXCK.L*2.-GV*Ll<K+2, 

2 lh ) ABS(GVALl(V,L) UGToRMAX) GO TO 10*1 
IF( ABS (GVALl (K » L ) ) o LToRMlN ) GO tO 1 1 

GO TO 149 „ 

1001 GVAL 1 ( K » L ) =Co 0 ' 

148 lPiE:§T:Kl! 88 J8?5 

I F ( GVALl f K t L + l 1 • LTo'^o 1 ) GO TO 75 

I F (GVALl ( K-l » L + l ) « LToOo 1 ) GO TO 75 
CALL MOLE ( MM ? K ? L ? X SC ?YSC? DIV ) 

^tilH:!:US!IStl'(5:l:l"^xSiXvM.uK-i,L*i»-Gv*Li(K + i, 

88 to iool 

GO TO 149 

1002 GVAL1( K»L)=CoO 

75 IF(LoLTo3) GO T 0 149 

82 I F ( Ko GT o N2 ) GO TO 185 ■ 

IF(GVAL 1 (K+2? L ) oLToOol ) GO TO *85 

165 IF(KoEQ.NMR) G0T085 _ „* 

I F ( GVAL 1 ( K+l- L -1 ) o LT 0 9. 1 ) GO TO 85 

N=N+2 
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85 IFtGVALl(K,L-2).LTo0ol) GO TO 149 

- t ° 91 

GVALI^k'i tl=GV&Ll(KTL) - "^AT?0^(GVALl(K + liL )-GVALl(KiL-l 
2 ))/(2oO*XR) dm i m ) GO TO 1007 

l?JftiUgv&tllK:tn:bT:5H£xi go to i<*7 

GO TO 149 

1C07 GVALUKtL =CoO 

91 lFIGViLl IK-1?L-1 NLT.Co 1 I GO TO 1« 
l lP( 1 GV»LnK,L-l).LT = -).ll GO TO 96 

94 1fIn.LT. 14) GO TO 96 

l^AtliK’h) =lviLi <I;l)"-Ia 1-?GVAL1<K+1,U-1 1-GVALKK-l. 

b ts mi 

GO TO 149 

1006 GVAL1 ( K » L) ■ « ’ it 0 o l) GO TO 149 

96 TC 199 

|X5Hi?:b\=§«ALl ti:L) _ -RATl5*tGVAH<K,l.-2t-GVALl»K-2,l. 

l§ *8 mi 

GO TO 149 

1005 GVAL1 ( K, L) =CoO 

140 IF(WI«EQ®1)G0 TO 142 G g tq ^49 

5 G^Li ( o|gV A?1 IS? L > -GV AL1 ( K , L + 1 ) ) / 

2(YSc£(GVALl(K,L-l)-bVAll(K,L+l) ) > 

2 £^\ Jl f 1 *GV AL 1 « K-2 . U 1 «VAL1 c K . 

2GVAL1 ( K,L-1 ) +GVAL1 <K t L 2 ) J / 

IF( ABS( ROX oGTog-. . ) GO TO 

IF ( ABS ( ROX) o GE o RAT I 0) G J iu 

RA X=ROX 
GO TO 1022 
1021 RAX=RAT I 0 
RAY=Oo r ' 

1C22 fp!Kl!§S8:§tGo*?ri?-2i 

ray=-c*b 
GO TO^ 1024 

1G 24 I F ( Mo LT o RJ1 ) GO TO 157 ^53 

153 g8 T8 ic .54 

p=ray/rax 
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1C53 I F ( ABS ( RAX) o LTo EPS I ) GO TO x05j> 

P=EPS I / RAX 
GO TO 154 
10 54 P=RAY/EPSI 
GO TO 154 

1 1 54 WRI TE ( 6 , 323 ) RA X , RA Y , R ATI 0, K , l 

1 57 I F ( GV A L 1 {*K » L ) ?GT I RM I nT GO TO 150 
GVAL 1 ( K. , L ) =0o 0 

150 IF( ]T JSCoEQoMESH) go to 170 
N= 1 

SUM=OoC 

DO 162 L=1,MMR 

YG=YSC*L 

DO 162^K=1,NMR 

IfTgVAL 1 (KtL)«LToO«l) GO TO 162 
0BS1 (NR+N) =GVAL1 (K, L) 

YO (NR+N ) =YG 
XC ( NR+N ) =XG 
N=N+1 

SUM=SUM+GVAL1 ( K , L) 

162 CONTINUE T , kt , . . 

AVER=SUM/ (FLOAT ( N-l ) ) 

NR=NR+N 

JSC=JSC+1 

NMR=2#NMR 

MMR=2*MMR 

XSC=XSC/2«C 

YSC=YSC /2» 3 

WRIT E ( 6 5 2 S3 ; X SC • YSC , NMR t MMR , NR 
GO TO 2 

170 DO 173 L=1,MMR 
1 3 = MMR + 1 
I =I3-L 
YG=YSC~ , 'L 
D0_173 K=1 , NMR 

IfTgVAL 1 (KtL)oGToCol) GO TO 172 

172 XGA=XG*£osTthIt)+YG*SIN(THET) 

1 YGA=-XG^SIN(THET)+YG*COS(THET) 

XG A=R MI NG-XGA 
YGA=RMI NL+YGA 
G(K, I ) =GVAL1 ( K » L ) 

173 CONTINUE 

WR ITE ( 6 , 255 ) 

WRITE(6t259) ( ( GVAL 1 ( K , L ) , K, L , K=1 , NMR ) , L= 1 , MMR ) 

IDR= INT-1 

DO 171 I =1 , I DR 

CL(1)=SINT 

171 CL( I + 1 )=CL( I ) + DELTA 
LTG(l)=IS«E0o2 
LTG(2)=JS.EQo2 

CALL 3 C0NTUR C (GVAL1,NW,MMR,5C,CL,1NT,T1TL6,9,21,LTGI 

111 VALUES FOR GRID POINTS,/! 

256 FORMAT (' <’'* , 30X,4HTEMP,5X,1HK,5X,1HL) 

258 FORMAT ( 2X, F15o2,3I5 ) 

P; PglHtT 1 ! ^ t3dI^TX I :IS|.§NgT i SdI : 8 Jc, 5 HDEPTH , 5X 

2 4 HT E MP t 2 X 1 1 1 H S CUN D VELOC,/) 

272 FORMAT ( I5,Fl‘,o 2,1 5,5Fl-jo^) 

273 FORMAT ( 12X,3F15©4) . 

274 FORMAT ( 20X, I 5 , F 1 0 « 2 , 1 5 , 5F 1 w ® 2 ) 
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280 FORMAT* 2F1C. 0.315) 



FORMAT ( * R' T 20X * 36HPCS I T ION OF THE GRID REFERENCE POINT 
Ifl FORMAT < 26X»8HLATlTUDEtSX » 9HLONG 1 TUDE j / I 
292 FORMAT ( I5fFiGoO,I5t F l'oC ) 

III 



4X » 1 



321 

322 

323 

324 



FORMAT ( ' i • I^oXt29H RELATIVE POSITION liN 

FORMAT ( 20X » 1HX » 15 X » 1H Y » 10X » 9HP ARAMETE R» / ) 

FORMAT ( lf'X t 3F1 5. 3, 2 I 5 ) 

325 FORMAT! 5ll?2GX*lHHESH,2X,7HHINIMUM,3X,7HMAXlMUM,2X,8HM 

\\% FORMAT^”:Io^32HiNTERPOLATED VALUES FOP THE GRID,/) 
350 FORMAT(12X,3F15.5) 

Ki^^X 8^T 0N,IX,5HPARAM,1X AHGRIO.2X.3HIN 

IS* • X,7HEPSIL0 

i FORMAT ( 2X, 315, F I 5. 2, 5X , I 5, El . 



_ _ . 2 ) 
POINT 



HAVE RECIVED AN AVERA 



4 30 . — 

||? FORMAT ( !xi 4RHTHE 'FOLLONI NG 

898 F 0 R M AT ( 6 A 8 ) 

899 FORMAT (315) 

STOP 

END 

c PROG RAM^F OR ^ lI A ST^ SQUARE^S UR FACE INTERPOL AT. ON UP TO SECON 
COMMON/ PAR AM/X ( 20 ) » Y( 20 ) » Z ( 20 ) 

DIMENSION C ( 6 » 6 ) ,B(6) 

I F ( No LTo 3 ) GO TO 2^ 

M= 6 

DO 2 1=1.6 
R f T ) = 0*0 
Du 2 K= 1 1 6 
2 C( K» I ) =OoO 
C( 1 , 1 ) =N 
DO 5 J=1 » N 
C(1,2)=C(1,2)+XJ J) 

C( 1,3) =C ( 1 ,3)+ v ( J) 

C lr4)=C(l,4)+X( J)**2 
C ( 1 ,5) =C ( 1 t 5) +X ( J )T-Y( J ) 

C 1 fc)=C(i!6)+Y( J)**£ 

C( 2 ,4) =C ( 2 ,4) +X( J ) **3 
C(2 5)=C(2,5)+Y( J)*X( J )**2 
C(2,6)=C(2,6)+X(J )*Y( J )* Jp 2 

C ( 4 1 i ) = C ( 4 ! 4 ) + X ( J U * A 

EifctUS! 3:l!n«i)2i2il(?(j)**2) 

II! 

B ( 1)=B( 1 )+Z(J ) 

B ( 2 )=B ( 2 )+Z( J )*X( J ) 

B( 3)=B( 3 )+Z ( J )*Y ( J ) 

B4=B(41 + Z(Jl*X|J)f2 
B ( 5 ) = B ( 5 ) + Z ( J ) * X ( J ) *Y( J t 
5 B(6)=B(6)+Z(J)*Y(J)**2 

C( 5 , 5 ) =C ( 4 , 6 ) 

C ( 3 » 5 > = C ( 2 , 6 ) 

C(3,4)=C(2,5) 

C ( 3 , 3 ) =C ( 1,6) 

C(2,2)=C(1 ,4) 

C ( 2 , 3 ) =C ( 1 , 5 ) 

DO 7 1=1,6 
DO 7 K= 1,6 
7 C(K,I)=C(I ,K) 

I F ( No GT • 5 ) GO TO l v 
M = 4 

DO 9 1=1,6 
C( 4, I )=C( 5, I ) 
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9 C ( 1 1 4 ) =C ( I » 5 ) 

10 CALL SIMO (C»B»M t KS) 

POLY = B ( 1 ) 



2UX(l)ix(2))i(Y(l)-Y<2t)) 

22 CHiu(M-Z(2nHtn/"(i)-Y°2W ' 4 

GO TO 28 7 , /? 0 

24 POLY= ( Z( 1)+ Z(2))/2oC 

?h cM=(9(t J-Z12) )*Xm/(X(l)-X<2)) 

28 XY= X lU*2 + Y(l5**2)**0.5 

CS=lo57C3*XY/XR 
POLY=Z( 1 )-CM 
42 RETURN 
END 

SUBROUTINE MOLE (MM ,K , L t XSC , YSC » D I V ) 

COMMON/ GRID/U ( 50 t 5 ... ) 

FAC1 = ( X SC/ Y SC **2 

5SS?^l52”S4KI)*f.* JH 

1 m'?K T L> -Z’- ? i 3(K + l ’.Lil? - uTk;'l + 2 l-U( K + 2.U) -UlK+2,1+2) 

1 DIV=^ulK+lfL+l )-U(K+2tL>2) 1 / X^ 

2 M*! kIl ) = -FAC2*(U(K + l»L + l) -2o^*U (K)L+l>+U(K-ltl-+l) J + 2o C* 

**2 U(K »L + 1 )“U( K» L+2) /v< - r 

dVv=?U(K,L+1 )-U(K»L+2) )/YSC 

, fjPJp.L;! ■0*U(K-l.l+D-UtK,L+2l-U(K-2.LI-U(K-2.l.+2) 

3 DIvi(U( K-l , L + l ) -ill K-2. L + 2 U /XR 

„ 8?Kl2)15MCl.(U«K-l,L*l)-2.0*U(K-l.L>+U(K-l,L-nl*2.Cv 

2 OIV-7u(’k'ia)-U(K-2,L))/XSC 

s S?KT^i2.0W(K-l.L-l»-y.K T 2.U-U,K.L-2l-U«K-2.L-2» 

DIV=(U(K-l,L-l)-U(K-2,L-2))/XR _ 

6 §?KT2)i2«C2.( U <K + l.L-l)-*.0*U«K.L-l> + U<K-l,t-lM* 

2 U(K.L-n 7 UtK,L-2l 

DI V=( U( K ,L-1 ) “U ( K, L-2 ) ) / » o L 

7 U?k!l)=4oO^U(K + 1,L-I)-U(K+2,L)-U(K,L-2)-U(K+2,L-2) 

DIV=(U( K+l »L-1 ) -U( K + 2* L-2) ) /XR 

8 S?KT?li2F»Cl.(U«K*l,l«»-2.e.U(K t l,Ll + UtK t l,L-in*2.9- 
2 oiv=tu(k+l'a)7u(k+2,Ln /xsc 

9 S?KlE. 1 4u(K*l.U +F ACl*(U(K,L + l) + U(KA-in^(K-l,LH/( 

2 2o^FACl+2oC) 

10 RETURN 
END 
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